Supersymmetric field theory of local light diffusion in semi-infinite media 



o 
o 

(N 
o 

Q 

(N 



I 

CZ3 



I 

o 



oo 
m 

O 



Chushun Tian 

Institut fiir Theoretische Physik, Universitdt zu Koln, Koln, 50937, Germany 

(Dated: February 2, 2008) 

A supersymmetric field theory of light difi'usion in semi-infinite disordered media is presented. 
With the help of this technique we justify-at the perturbative level-the local light diffusion proposed 
by Tiggelen, Lagendijk, and Wiersma [Phys. Rev. Lett. 84, 4333 (2000)], and show that the 
coherent backscattering line shape of medium bar displays a crossover from two-dimensional weak 
to quasi-one-dimensional strong localization. 
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I. INTRODUCTION 

The Anderson localization of light has been one of the 
most fascinating phenomena in condensed matter physics 
since the mid-eighties^'^. Like electron systems this phe- 
nomenon finds its origin in coherent multiple scattering 
which slows down diffusion of photons and eventually 
brings them to stop. Parallel to studies of disordered 
conductors the subject in this field ranges from light lo- 
calization near or far below the mobility edge in bulk (in- 
finite) systemsii^ (where low-energy photon motion en- 
joys the translational symmetry) to their detection such 
as transmission measurements in the slab geometry (e.g., 
Refs. 



A unique subject of localization in optical (and 
other classical wave) systems is the enhanced coherent 
backscattering (CBS) phenomenonii^. In this subject 
the issue of semi-infinite geometry is heavily addressed 
because the CBS line shape is responsible for by optical 
paths near the vacuum-medium interface. Although it is 
well known that in the weak disorder region, i.e., I ^ A 
{I the mean free path and A the wavelength) incident 
photons enjoy diffusion as in bulk media^'^, in the strong 
disorder region / < A the role played by the leakage at the 
interface has been of long term interests^ and, still, re- 
mains in the central position of CBS studies, particularly 
to forecast or observe the CBS line shapei^iiiii^ii^. 

Pressingly, in the latter region strong localization 
emerges in the bulk and a new scale namely the lo- 
calization length ^ appears. On the experimental side, 
there has been increasing evidence indicating that in- 
side the boundary layer of thickness ^ ^ the photon 
leakage at the vacuum-medium interface strongly inter- 
plays with strong localization^i^iiiii^. On the theoret- 
ical side, some time ago exact solution of semi-infinite 
one-dimensional disordered chains shed the light on the 
existence of so-called radiative localization states in the 
boundary layer—, which lead to anomalous slowing decay 
of reflected (backscattered) incident light pulse a^^i^^ . 

Recently, in an insightful theoretical worki^ it was re- 
alized that (in three-dimensional disordered media) in- 
side the boundary layer the translational symmetry of 
low-energy (hydrodynamic) photon motion is strongly 
destroyed resulting in the so-called "local diffusion" . 



Remarkably, constructive wave interference renders the 
static diffusion coefficient depending on the distance 
from the interface. Consequently, the local diffusion 
was found to lead to a rounded CBS line shape resem- 
bling that observed experimentallyii and, thus, might 
overcome the conceptual difficulty of earlier theoreti- 
cal proposal^. Surprisingly, the dynamic generalization 
of the local diffusion equatio n^^i^^ provides an expla- 
nation of some key phenomena observed in quasi-one- 
dimensional microwave experiment*^, and well captures 
anomalous slowing decay of reflected incident pluses in 
both quasi-one-dimensionall-^'^^'^^ and three-dimensional 
disordered mediai^. Moreover, such novel prediction-the 
position dependence of diffusion coefficient-seems to have 
been confirmed by a very recent experiment^"*. 

Despite of these progress theoretical investigations 
on local diffusion are, however, restricted on the self- 
consistent diagrammatical metho d^^'^^i^^ . Thus, an in- 
tellectual challenge is to seek for the genuine micro- 
scopic origin underlying this novel concept. This is, in- 
deed, the purpose of this work. The last few decades 
have witnessed spectacular success of applications of su- 
persymmetric field theory to various disordered systems 
in the absence of interactions (Such condition is per- 
fectly satisfied by optical systems.)^". Among them 
there are few exact nonperturbative results for quasi-one- 
dimensional disordered wires such as density-density cor- 
relation function (in the infinite geometry )2i and trans- 
mission statistics—. They allow one to make important 
insight on the strong localization. Most importantly, for 
periodic disordered media by using the supersymmetric 
field-theoretic method a local light diffusion equation, 
similar to that proposed in Ref. Il2l . recently has been de- 
rived at the microscopic level^^. In view of these it is nat- 
ural and inevitably necessary to proceed along the same 
line to explore the concept of local light diffusion and its 
effects for more general-fuUy disordered-media, which 
differ drastically from the former~ from both physical 
and technical view. 

The main results of this paper are as follows, (i) We 
present a field-theoretic proof showing that, in contrary 
to the conjecture of Ref. 19|, no scaling behavior exists 
inside a layer of thickness ~ I extrapolating into the vac- 
uum, (ii) We justify-at the perturbative level-the local 
diffusion equation proposed in Refs. [i^ITtI . (iii) We an- 
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alyze signatures of the static local diffusion in the CBS 
line shape. It should be stressed that in this paper the 
suspersymmetric field theory is treated perturbatively, 
and the nonperturbative treatise will be reported in the 
forthcoming paper. 

The rest of this paper is organized as follows. In 
the next section we produce nonlinear supersymmetric 
a model in the context of optical systems. Most impor- 
tantly, we derive the boundary constraint satisfied by the 
supersymmetric matrix field. The supersymmetric field 
theory is then applied to the two-dimensional medium 
bar. Sec.|TTT]is devoted to exploring states residing deeply 
inside the semi-infinite medium bar (namely far away 
from the interface) by investigating renormalization ef- 
fects of infinite medium bar. In Sec. IIVI weak localization 
in the semi-infinite medium bar is studied, where the 
general dynamic local diffusion equation is justified. The 
static limit of the local diffusion equation is studied in 
Sec. |Vl In particular, the weak localization correction to 
the bare diffusion constant is explicitly calculated, and 
its signatures in the CBS line shape are analyzed. We 
conclude in Sec. |Vl] and give some technical details in 
Appendix [XllEl 



II. SUPERSYMMETRIC FIELD-THEORETIC 
FORMALISM 

In this section a supersymmetric field-theoretic formal- 
ism is presented for light scattering in semi-infinite dis- 
ordered medium. 



A. Nonlinear a model 

We first show that as interactionless electron systems 
low-energy photon motion in bulk disordered media is 
well described by the nonlinear a model. The derivation 
is rather standard^S. Here we only outline the scheme 
with an emphasis on the main difference, while refer the 
reader to Ref. 20. for the details. 

In the present work for simplicity the scalar wave will 
be considered. The wave propagation in a bulk disor- 
dered medium is described by the Helmholtz equation: 



{V^ + ^f [l + e(r)]}i?(r) 



(1) 



where the field E has the radiation frequency 17 (velocity 
c set to be unity), and j(r) is the source. Here the fluctu- 
ating dielectric field e(r) has zero mean and is distributed 
according to the Gaussian ^-correlated law: 



(e(r)e(r')) - A5(r-r'). 



(2) 



The Helmholtz equation resembles the Schrodinger 
equation with the Hamiltonian now read out as H — 
— — V.^ e(r) . As usual we may introduce the re- 



tarded/advanced Green function 602^ defined as 



{n 



S{r- 



(3) 



where r2± = ± iO~^ . The electric field and the source 
are related via E{r) = J dr' 6*^2 (r, r')j(i"') ■ We may 
further introduce the diffuson and the cooperon y'~^ 
propagator defined as 



y^{r,r';u;) ^ Gf^+^+/2)2 (r, r') Gf^_+/2)2 (r', r) . 



:yC(r,r';c.) ee Gf, 



(r.rOGf^_.+ /2i.(r,r') (4) 



/2)2 



with uj^ — uj + iO^ and w <C li , where the overline stands 
for the average over random dielectric field. These two 
propagators describe elegantly the light propagation over 
large scales. 

The propagators above are represented in terms of su- 
perintegrals. For this purpose we define a supervector 
field -0 : 




,m = 1,2 



(5) 



with S"s (x's) complex commuting (anticommunting) 
variables, where the superscript 1 (2) refers to retarded 
(advanced) Green function, and its charge conjugate 
ijj = ■0^A . Here A is an 8 x 8 supermatrix: 



A 



1 

-1 



|bf , 



(6) 



Hereafter supermatrices are defined on the re- 
tarded/advanced (ar), bosonic/fermonic (bf) and 
time-reversal (tr) sector. Then 

:y°(r,r';c^) 



-4 J d[V']Vi(r)V;i(r')V^(r')^^(r)e-^[V'.^l, 
3^^(r,r';^) 

-4 1 d[^]V'i(r)V:;i(r')0^(r)^^(r')^^^[^. (7) 
i I tpir) [-Ho - n^e{r) - n uj+A] V'(r)dr . (8) 



Here 



with Ho = — — , where the term is omitted since 
Lu ^ . Performing the average we arrive at 



(9) 



The quartic term is decoupled by the standard 
Hubbard- Stratonovich transformation. Introduce an 8 x 8 
supermatrix field Q{r) conjugate to ^^^^^2) V'('') ® V'(r) • 
Here iV(ri^) is related to the photon density of states per 
unit volume z/(0) by viVl) = 2nN{9?) . Then, 



exp 
= /cxp 



7riV(f]2) 



Q^]dr 



(10) 



D[Q]. 
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Substituting it into Eqs. ^ and ([9]) and integrating out 
the tp fields using the Wick theorem, we obtain: 



(11) 



(str[fc(l + A)(l - r3)Q(r)(l - A)(l T r3)kQ{r')]) 



with 



1 

-1 



bf 



(12) 



and Tfc the Pauh matrices defined on the time-reversal 
sector. In Eq. pT|) the following average is introduced: 



(P[Q]) = J D[Q] P[Q] e 



-FIQ] 



where the action F[Q] is 



AQ2 



(13) 



(14) 



F[Q] = J drstr 

i In [-iHo - inuj+A + TTAN{n^)Q{r)] | . 
Minimizing F[Q] gives the saddle point equation: 



Q 



1 



1 



{-iHa - iiluj+A + TrAN{n^)Qy 



Go 



(15) 



In the limit O » uj ,TTAN{n'^)/n , Eq. ^ gives the 
saddle point as Q{r) = A. 

So far the derivation is exact. Fluctuations analysis 
around the saddle point may be performed for Eqs. (jl3p 
and p4p . Yet, we could not proceed further and on ly g ive 
the results here, instead refer the reader to Ref. for 
all the details. First, after standard procedure the mean 
field approximation namely Q{r) — A gives the averaged 
retarded/ advanced Green function as 



G^^^(r,r') = {r\{n^ + ±nrANin^)}~^\r') . (16) 

The imaginary part of the self-energy gives the elastic 
mean free path which is 

" nAN{n^) ' ^^'^^ 

and has the Rayleigh form, i.e., / ~ il^^'^+i) . 

Then, with 0/^1 taken into account the action is 
simplified to be F[Q] = J dr C[Q] , where (From now on 
we set = ^{Q) to shorten formula.) 

^[Q] - ^ str [Do [dQf + 2iu;+AQ] (18) 
o 

with the bare diffusion constant Dq — l/d. Here 
Q{r) — T(r)AT~^(r) describes Goldstone modes with 
T(r) a matrix field taking the value in the coset space 
C/(2, 2/4) /C/(2/2) X ?7(2/2) reflecting the orthogonal sym- 
metry. An explicit parametrization of T will be given in 
the next section. 



B. Q-field constraint at the vacuum- medium 
interface 



The action F[Q] obtained above is invariant under the 
translational symmetry, which is broken in the presence 
of the vacuum-medium interface. The broken transla- 
tional symmetry may profoundly affect light propaga- 
tion. Experience in mesoscopic physics shows that to 
take into account the vacuum-medium interface effect one 
may impose some appropriate boundary condition on the 
Q-field in the field-theoretic formalism. However this is 
a nontrivial task and in meseoscopic physics investiga- 
tions so far have been restricted on interface structures of 
quasi-one-dimensional disordered wires and small quan- 
tum dot o^°i^^'^^i^^ . In this part we switch to optical sys- 
tems and study the vacuum-medium coupling where the 
interface may be infinite and bear arbitrary geometry. 



1. The vacuum-medium coupling action 

Though the derivation below may be generalized to 
arbitrary dimension to simplify discussions we will focus 
on the two-dimensional case. Let us suppose an arbitrary 
curve C which divides the space into two disconnected 
subspaces V_ and V+ , i.e., = v_ U V+ U C and V_ n 
V+ = . We are interested in light propagation in some 
subspace say V+ described by an effective Green function 
5j^2^(r , r') , which is identical to G^2^{r , r') for r , r' e 
V+ . To study such Green functions for r , r' € V_ we 
introduce auxiliary Green functions gQi^iic , r') satisfying 

{nl-H}g^i^{ry) = 5{r-r'), (19) 

Then our starting point is the following theorem due 
to Zirnbauer— and refined by Efeto\*2£, which was 
originally established for description of coupling between 
leads and mesoscopic devices. The theorem is stated 
as follows: (For the self-contained purpose the proof 
tailored to the present context is given in Appendix [Xl ) 

For r , r' G V+ the Green function Gqi^Itc , r') solves 

{nl-H± iB} g^i^ir, r') = <5(r - r') , (20) 

r e C, 



9n(r)^si2'^(r,r') 



rec 



0, 



where the normal unit vector n(r) at r points to V+ . 
Here 

(B/)(r) = / dr'Im[S(r,r')]/(r'), 



S(r,r') = 9„(r)9„(r05fl2(r,r'),/orr,r'eC. (21) 

The effective Hamiltonian for the retarded (advanced) 
Green function is H ^ iB . Remarkably, it is non- 
hcrmitian due to the escape from V+ into V_ through 
C. 
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In the present case the vacuum-medium interface 
namely the curve C is a straight Hne. To proceed we 
choose the coordinate system {r± , z) with the z {r±)- 
direction perpendicular (parallel) to the vacuum- medium 
interface. The vacuum fills the space z < where no di- 
electric scatterers are available. For technical reasons we 
assume that the dielectric scatterers located at (r^ , 0) 
{r]_ < < • • • ) are uniformly (in the statistical sense) 
distributed with the distance between nearest scatterers 
li = r\ — r^^^ order of L C is located at z = 0~ and V+ 
(V_) is set to be the medium (vacuum). 

Taking into account the boundary condition specified 
in Eq. (fT9)l we find that the Green function g^a is 



(22) 



dk- 



sin(/cz) sin(fcz') 



dk I 



Upon the substitution of Eq. ([2211 into Eq. ([2T|) the op- 
erator B is simplified to be 



dk I 



(23) 



\k±\<n 



^ / dr'^ ^n^-kl 6^'=^ ^"^""^^/(rl). 



Repeating the derivation of Sec. Ill Al with the effective 
Hamiltonian H ^ iB we arrive again at Eq. except 
that the action is modified according to F[Q] ~^ F[Q] + 
-Fintor[Q] with 



FintortO] - strln {l - BAGo] 



(24) 



Expanding the logarithm and substituting Eq. (1231) into 
it we obtain 



1 



■ str ■ 



oo 

E 

.n=l 



BAGo 



where the supertrace str includes the integration over r± , 
and fJo(r,r';Q) exponentially decays for |r — r'\ ^ I ac- 
cording to Eq. (I15|) . To calculate Eq. (|25|1 we introduce, 
for arbitrary i , the auxiliary variable = — , 7'^ G 
[r^"\ r\_] . Following Ref. HO in the layer < z < I 
the Green function C/o(r,r';(5) = Goir^ — Xi,z,r\ — 
x[, , z', ; Q) may be approximated by 



Go{f 



Xi, z, rj_ , z , , Q) — 



2(5,.' 



N>1 



(26) 



''^^(1 



"PIn/I, i^i)'PlN/U (^^) COs(fcz) COs(fczO 

') - fc2 - fc2 + iT:l\N{9?)Qi ' 



where the longitudinal wave function is determined by 
the boundary condition of Eq. ([20]) . Here and 1 -I- 
stand for the Q- and the dielectric field, respectively in 
the regime: [r^^ , r^] x [0 , . They are considered to 
be a constant (matrix) since both Q and e varies over 



the scale I . Moreover, the transverse component Lp\^ is 
defined as 




(27) 



Substituting Eq. ^ into Eq. with the help of 

the following identity: 



TtN\ /, TTiV 



(28) 



we obtain: 

Fintcr[Q] = 51 E Strln + (29) 



0<Tf <i 



with 



(30) 



Taking the advantage of large channel number 
^o<*^<i ^ f^^j/TT we may simplify it to be (The de- 



tails are given in Appendix [BP 



4e 



To(z) str (AQ, 



(31) 



by assuming that a\,^ does not depend on k±_ , i.e., a\,^ 
a' . Here 



4a* 



< 1 



(32) 



(25) is the well known transmission coefficient of electromag- 
netic wave^l. Passing to the continuum limit: U 
J dr± we rewrite the action i^intor [Q] as 

^^intcr[0] = J dr±_ Toir^) str [AQ(rx, z = 0)] . (33) 

In the quasi-one-diniensional geometry the summa- 
tion over i is suppressed, and the coupling action 
-Flntcr[(3] namely Eq. ([3T|) recovers the one obtained 
previously^"'-^^. For d > 2 although to generalize the 
derivation above is straightforward, the coupling action 
may be obtained by simple physical arguments below. 
Notice that the coefficient of Eq. (pij) allows a simple 
physical explanation^: According to Eq. (fTB|) the (sin- 
gle) photon Green function decays over the scale / . Sup- 
pose that the medium is partitioned into boxes of volume 
l"^ , then the states (denoted as /i) in different boxes are 
uncorrelated. The box states neighboring to the inter- 
face may be translated into the vacuum state (denoted 
as a)-so-called lead channel in the terminology of meso- 
scopic physics. The couphng strength is oc J2f_i Wa^W^a 
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with Wafj. the scattering matrix element, which scales as 
l'^~^/A with A the interface area. Thus, although the to- 
tal channel number is cx A^'^^^ , the number of channels 
to which the interface box state is transmitted (denoted 
as Nd) is much smaller Nd ~ Afl'^-^ x I'^-^/A = {QlY'^ . 
More precisely, Nd may be found to be 



Nd = 



(2n 



dk I k 



d-2 



(34) 



For d = 2 this gives N2 — VII/tt namely the coefficient 
of Eq. (pij) . For arbitrary d with the replacement of 
N2 Nd (and I — > h) in Eq. (|3ip the vacuum-medium 
action becomes 



-Winter [Q] 



2^ Z Tq[i) str (AQi) 



4 ^ TT 

■i 



dr,5cro(r) str[g(r)A] (35) 



with 



N, 



Jd-l)/2 



{2n 



(36) 



In the last equality of Eq. p5|) we again pass to the 
continuum limit, and the operator Sc is defined as 
Jdr6cf{r) = Jdr^f{r^,z = 0). 



2. Boundary condition 

We then come to derive the boundary condition sat- 
isfied by Q . For this purpose we employ the so-called 
boundary Ward identity^^. It states that an arbitrary 
local observable say P(r) (with r inside the medium), 
which is expressed in terms of the average of the func- 
tional VlQir)] namely 

P(r)= /7:>[g]7'[Q(r)]e~^^eR+'^"-£[Q(r)]-^F.„,e,[Q] ^ (37^ 



must be invariant under an infinitesimal boundary rota- 
tion below: 



Q ^ e-''Qe''^Q-[R,Q], 



(38) 



R = 





n{r±,z = 0) 







Notice that the boundary rotation alters neither 
^[(^(r)] nor £[Q(r)] for r inside the medium. The bound- 
ary Ward identity then demands Ji^intcr = , i.e.. 



■To[Q,A] 



0. 



As 7?. , 72. are arbitrary this requires 



(!Qd,Q + To[Q,A] 







(40) 



to be met, where the subscript _L stands for the offdiago- 
nal component in the retarded/advanced sector (thereby 
anticommunting with A) and / = 2'kvDq/ [Nd^'^^^) ■ 

Eq. (HDl) is the field-theoretic version of the radiative 
boundary condition— As we will show in Sec. IIVI 
it describes that the low-energy coherent dynamics pen- 
etrates into the vacuum with the extrapolation length 



^^2To ^Twd^Yo 



(41) 



Notice that it is proportional to the inverse transmis- 
sion coefficient in agreement with Ref. [29l For d = 3 
in the case of perfect transmission, i.e., Tq — 1 the ex- 
trapolation length is C = I Hn agreement with Ref. 
and is closed to the one that obtained by solving Milne 
equation^ which gives C = . Traditionally the radia- 
tive boundary condition is imposed to diffusion equation 
to mimics the leakage at the interface^^'^''''^"''^^ and is 
justified for one-dimensional discrete random walk^. 



III. TWO-DIMENSIONAL 
RENORMALIZATION EFFECTS OF INFINITE 
MEDIUM BAR 

In the rest of this paper we will apply the super- 
symmetric field-theoretic formalism to the semi-infinite 
two-dimensional medium bar (with the width a 3> l)^ 
where in the transverse (/?-) direction the photon mo- 
tion is confined. The purpose of this section is two- 
fold: On the physical side, we wish to explore how a 
finite width affects localization in the bulk which dif- 
fers inessentially from localization in an infinite bar. Ac- 
cordingly, through out this section the action reads out 
as F[Q] = J^^dz J^dpLlQ]. On the technical side, 
by presenting some details we wish to address the dif- 
ference of calculations between semi-infinite and infinite 
bar, which originates at the fact that in the former system 
the translational symmetry of low-energy modes, i.e., the 
Q-field is broken. 

Following the standard strategy we factorize the T- 
field into the slow and fast mode in terms of T = Ty.T^ , 
where rotations T> (T<) involve spatial fluctuations on 
short (large) scales. Substituting it into the action we 
then obtain: 



00 fa 



00 Jo 



dz dp Sir {Da {{dQ>)^ + 



(42) 



4g>9Q> + Q>]2) + 2^cJ+g>T-lA^<}, 

where Q> = T> Ar> ^ and $ = T^^dT^ . Integrating out 
(5> results in an effective action of (5< . 
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A. Parametrization of fast modes 

To work out the strategy outlined above we set r> — 
1 + iWy with parameterized by 



B> 

By 



(43) 



Since photons are confined the current vanishes at p — 
0,a, i.e., dpWy{r)\p=o — 0. We may thus introduce 
the Fourier transformation 



J\k.\^Ko n>no 



(44) 



where kg and imo/a are ultraviolet cut-off of longitudinal 
and transverse wave number, respectively. In Eq. (j43[) the 



matrix By has the structure as 



By = 



a la 
rj ih 



bf 



(45) 



with 



(Tl (72 



hi 62 
b*2 hi 

V2 V*i 



(46) 



where a's, 6's (cr's, ?7's) are complex bosonic (Grassmann) 
numbers, and the charge conjugation transformation of 
a matrix M is defined as 



M^CoM'C^\ Co = 



-tT2 
Tl 



(47) 



Straightforward calculations justify the useful identity: 
M1M2 = M2M1 . 

Importantly, Wy satisfies the following relation: 



Wy = KWlK , K = 



1 
k 



(48) 



which, as shown in Appendix [Cl enforces the invariance 
of Wy under the charge conjugation, i.e., Wy — Wy . 



B. One-loop renormalization 

Now we study the one-loop renormalization. In doing 
so we expand Q> to quadratic order in Wy . Conse- 
quently the action separates into three contributions: 



F — Fs + Fp + FsF , 
where the slow mode action is 



(49) 



F5 - — / dz dpstT[Do {dQ<f + 2iuj+kQ<] (50) 

<5 J-oa Jo 

with Q< = T<AT< , the fast mode action is 



dz dpstiidWyY , (51) 
00 Ja 



and the slow-fast mode coupling is described by the ac- 
tion: 



/CO pa r 
dz dpstil Da {[Wy,dWy] ■ $ - i^AWyf - {<S>AfW^ - idWy ■ $ - i{<i>AfWy 
00 Jo I 



{iWyKT~^AT<_ + WlKT-^KT<) 



In order to calculate the general average with respect to the fast mode action: {■ ■ ■)f = /-D[VK>] (• 
employ the contraction rule below: 



2T:i'{Wyir)MWy{r'))F ^ -X>F(r,r')M, M = 
4:Tri'{Wyir)NWy{r'))F = X>F(r,r') (strTV - AstrATV) 



Af 12 ^ 


. ar 





) ' 


N = 


/iv" 


M1.2 


-iMfi 



(52) 



)e we 



r 



(53) 



with the fast mode propagator 
I?f(r,r') 
2 f dk ^ 

k\>ko „>„^ Do 



(54) 



and the Wick theorem. Notice that we assume that the 
fast mode propagator does not depend on the transverse 



cos ""(P-P') e'fe("-"') 
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center-of-mass coordinate, i.e., {p + p')/2 and the self- 
average over this variable has been performed. 

Performing the average of Fsf we obtain an effective 
action Fcff[(3<] = + {Fsf)f , where 



oo pa 



{Fsf)f = -^ I dzl dpstvidQ^)' (55) 

'-00 Jo 



with 



J f dk 1 



(56) 



Eqs. ([55| and ([56|l show that the one-loop renormaliza- 
tion results in the weak localization correction to the bare 
diffusion constant Dq : 



oc pa 



Fos[Q<] = — / dz dpstT{[Do + SD^^^]{dQ<)^ 

o J-oo Jo 



(57) 



with the two-dimensional weak localization correction 



C. Dimensional crossover of effective action 

In the high-frequency region, i.e., lo ^ D^^/a^ , the con- 
dition: koa/n ,no ~ ^ujo^/Dq ^ 1 is met. The photon 
motion is thus two-dimensional described by the action, 
Eq. ([57)1 . Furthermore, since 17/ ^ 1 the two-dimensional 
weak localization correction SD^"^^ is much smaller than 
Dq , the photon motion thereby is diffusive. 

In the low- frequency region, i.e., uj ^ Do /a? one may 
further enforce k^ — Q and no = 1 and thereby obtain a 
quasi-lD effective action of Q = Q<{z) which is homo- 
geneous in the transverse direction: 

/•oo 

Fcs[Q] = / dzCcs[Q] , 



C,s[Q] = —fitr[D,s{d,Qf + 2iu+KQ], (58) 



where the renormalized diffusion constant is 



D, 



off 



Dq 



n>l 



(59) 



Note that in the above the second term suffers logarith- 
mic divergence which, as usual, may be regularized by in- 
troducing the upper cut-off N which is order of ^ a/{nl) . 
As a result, 



D, 



Dq I 



I 



tt^vD, 



In- 







(60) 



Thus, in the low- frequency region: oj < Dq/o? the sys- 
tem is quasi-one-dimensional provided that the bar width 
satisfies a <^ le^ . For wider bar the system display 
two-dimensional strong localization which is beyond the 
scope of present perturbative analysis. 



WEAK LOCALIZATION IN SEMI-INFINITE 
TRANSPARENT MEDIUM BAR 



IV. 



The discussions of Sec. IIIII break down in the semi- 
infinite geometry due to the absence of the translational 
symmetry. In this section and the next we turn to study 
the vacuum-medium interface effect on wave interference. 



A. Simplified boundary condition 

In order to explore the physics implied by the boundary 
constraint Eq. ([40]) let us parameterize T in the same 
way as Eqs. (gS]), (gH]) and (gH]). (To distinct notations 
from those of Sec. IIIII we eliminate all the subscript > .) 
With the substitution of the parametrization and keeping 
Eq. (|40|) up to the first order in W we obtain: 



{Id, - 2Tq)B 

[Id, - 2Tq)B 







(61) 



implying B ,B ^ e^/'^ , z < with C = l/i^To) . Hence 
the low-energy Goldenstone modes penetrate into the 
vacuum of a depth ( then exponentially decays. That 
is, the optical paths underlying coherent multiple scat- 
tering do not cross the line located at z = —C, . 

From now on we assume that the interface is almost 
transparent namely Tq closed to 1 . In this case ( — . 
Since the mean free path I is much smaller than any other 
macroscopic scale we may safely assume that the crossing 
line where W vanishes coincides with C . Consequently 
the boundary constraint Eq. (|40)l is simplified as 



(62) 



QU=o = A. 



It is the action: F[Q] — J^^ dz J^dp£[Q] with the Q- 
field subject to this boundary constraint that we will use 
in the rest of this paper. 



B. Bare diffusive propagator 



Expanding Q in terms of W gives 
F[Q] = F2[W] + F4W] + 
where the Gaussian action: 



(63) 



and 



F2[W] = ^/ dz dpstT[DQ{dW)^ ~ iujW^] (64) 
2 Jq Jo 

/•oc /•a 

F4W] = ^ / dz dp {~2DQSti [{dWfW^] 
2 Jq Jq 

+iujstiW'^}. (65) 



From Eq. (|64[) immediately we obtain the same con- 
traction rules as Eq. (|53p except making the replacement: 



PF(r,r') ^I?(r,r';w) 



(66) 
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where the propagator solves the diffusion equation: 



r or r' G C 



= 0. 



(67) 



The boundary condition above is inherent from the con- 
straint Eq. (f62|) which imposes T4^(r)|rec = 0. 

Keeping the prefactor of Eq. ITTI) up to quadratic term 
we obtain the leading cooperon propagator: 



3^(6)(r,r';c.) 



(str{fc(l + A)(l-T3)M^(r) 
xk{l-K){l + T:,)W{v')})^^ 
(68) 



It is easy to see that the propagator above preserves 
the symmetry 3^^^^(r, r'; w) = 3^^*^^ (r', r; w) inherent from 



Eq. pTjl . Eq. is traditionally obtained by summing 
up all the ladder diagrams and imposing appropriate 
boundary condition^. 



Weak localization correction 



We proceed to calculate the one-loop correction to the 
bare propagator 3^^^^ . For this purpose we keep the W- 
expansion up to the quartic terms for both the prefactor 
and the action which, after straightforward calculations, 
gives the cooperon as 3^^ « 3^^^^ -I- dy'^ with 



<53'^(r,r';c.) = - 



7riV(r22 



n 2 



(str {fc(l + A)(l - T3)W^3(r)fc(l - A)(l + r3)W^(r')} 

+ str {fc(l + A)(l - T3)W^(r)fc(l - A)(l + T-i)W^{Y')] 
- str {fc(l + A)(l - n)W^{v)k{l - A)(l -I- T^)W^{v')] 

+ str {/c(l + A)(l - T3)W{v)k{l - A)(l + T3)W{v')} Fi[W]) . (69) 

First, it is easy to show that the third term in the right hand side of Eq. (pS)) vanishes. Second, as shown in Appendix [Dl 
the first two terms partly cancel the last term. Eventually Eq. (|69p is reduced into 



,53'^(r,r';c^) = - 



7riV(f72) 



n 2 



[-KvDo) / dzi/dpi(str{fc(l + A)(l-r3)Ty(r)fc(l- A)(l + r3)W^(r')} X 



str w {ri)W[vi)W {ri)W{vi) + {dW{vi) W{vi)Y + {dW{vi)Y W'^{ri)\ ) 

I 



F2 



(70) 



where the overbrace fixes the contraction and the deriva- 
tive acts only on the nearest W . We remark that Sy'^ 
vanishes when 3^^^^ is spatially homogeneous, which is a 
reflection of the flux conservation law or Ward identity at 
the one-loop level. Notice that Sy'^ vanishes if either r or 
r' belongs to the interface C . Such property is inherent 
from Eq. (fTTjl which vanishes upon sending either Q{y) 
or Q{y') to A. Using the contraction rules and integral 
by parts we further reduce Eq. ([70)) into 



on r°° ra 

Sy^'ivy-u) ^ j dz,J dp,V{r,,v,;cu) X 

dr,V{r,ri;Lo)d,,V{r',ri;Lo) (71) 



after tedious but straightforward calculations. Notice 
that ^3^*-^ preserves the symmetry: 6y'~^ {r , r' ; u) = 
Sy^ir',r;uj). 



D. Local diffusion equation 

Eq. dZIl) justifies that 3^^ = 3^^^^ + 53^'^ solves the 
following local diffusion equation: 

{-dD{r;uj)d-iuj}y'^{r,r';uj) ^ S{r~r'), 

y'^Uc - (72) 

at the one-loop level. Here D{r; u) = Dq + SD{r; uj) with 
the weak localization correction 

6D{r;uj)^~—V{r,r;u;). (73) 

The local diffusion equation differs from the traditional 
one in that the diffusion coefficient is position-dependent. 
It may be amounted to incompletely developed construc- 
tive interference between two counter-propagating opti- 
cal paths-which leads to the weak localization-near the 
boundary. Indeed, although deep inside the medium 
dD{r; lo) saturates recovering the bulk weak localization, 
at the interface it vanishes, i.e., 

dD{r;Lo)U=o = 0- (74) 
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Importantly, this is contrary to the theoretical proposal 
of Ref. which claims that wave interference democrat- 
ically renormalizes the diffusion constant appearing in 
both the diffusion equation in the bulk and the radiative 
boundary condition at the interface. 

Here several remarks are in order: (i) Higher order loop 
corrections preserve Eq. ([7^ . They affect the local diffu- 
sion equation by introducing higher order weak localiza- 
tion corrections which are also position-dependent. This 
peculiar property reflects the photon number conserva- 
tion law and is protected by Ward identity, (ii) In the 
presence of internal reflection namely To{r) (far) below 
I , (i) is no longer applicable because the simplification 
namely Eq. (|62|) breaks down due to large extrapolation 
length. In fact, Ref. falls into this case, (iii) The con- 
cept of local diffusion originally introduced in Ref. 12 at 
the static limit, i.e., w — s- together with its dynamic 
generalization^^'^* is now justified at the perturbative 
level. 



V. STATIC LIMIT OF LOCAL DIFFUSION 
EQUATION 

In this section we study the static limit: — > (For 
this reason below we suppress the argument uj in all the 
formulae.) of the local diffusion equation namely Eq. ([7^ 
for a bar with the width satisfying I <ti a <ti le^ " . 
In particular we will explicitly calculate the weak local- 
ization correction Eq. (|73p . and study its effects on the 
coherent backscattering phenomenon. 



A. Quasi-ID massive local diffusion equation 

In the static limit the weak localization correction 
Eq. l|75|) becomes self-averaged over the center-of-mass 
p and thereby is p-independent. That is, 



D 

5D{z)^ dpV{z,p,z,p) 

TTVa Jq 



(75) 



Substituting Eq. (j75l) into Eq. ([72]) we find that 
3^^(r',r) depends on p — p' , but not on the center-of- 
mass {p+p')/2 . Therefore, we may introduce the Fourier 
transform: 



y'^iz^z'^p- p') (76) 

i y^{z, z'\ + V ,') .OS 



n>l 



and insert it into Eq. ([7^ to obtain {q± = ^) 

{-d,D{z) + D{z)ql - iO+}y^^ {z, z') = 5{z - z') 



(^ = 0,z') = 0. 



(77) 



where 0"*" is infinitesimal positive constant. Eq. (j77p may 
be considered to be a quasi-one-dimensional local diffu- 
sion equation with a mass D{z)q'^ . 



B. Dimensional crossover of weak localization 

The weak localization correction Eq. (|75p then be- 
comes 



SD{z) = -^ 
nva 



Voiz,z) + 2j2 'D^iz,z) 



(78) 



Here Vq^^ {z, z') satisfies 



{Z?oHf + -z0+}l?gjz,z') - S{z-z'), 

Vq^{z^Q,z') = 0. (79) 



It is solved by (with the introduction of q±^ = \/q'^ — iO+) 



Vgjz,z') 



nva 



<1± 

1 e'-L'-e' 



-q±z 



,Z<Z' 



2? q± 

where ^ — niyaDo ■ Substituting it into Eq. (|78p gives 



5D{z) z a ^ 1 - e-^"^^/"^ 

^ ^ n>l 



Dn 



where the first term is the quasi-one-dimensional con- 
tribution, and the second term is the two-dimensional 
contribution with nixja standing for the transverse hy- 
drodynamic wave number. 

As expected at z = the weak localization correction 
(5-D vanishes. Away from the interface i.e., Z <^ z ^ a it 
may be approximated by 



bD{z) 
Da 



a z 
-In- 



(82) 



as shown in AppcndixlEl This suggests that in this region 
(even in the static limit) the two-dimensional low-energy 
motion dominates the weak localization. 

The two-dimensional contribution saturates at z ~ a : 



SDjz) 
Do 



z a X ^ I 



r!>l 



(83) 



where the second term is none but the bulk weak local- 
ization correction (see Eq. (|59p ) renormalizing the bare 
diffusion constant Dq . With this taken into account 
Eq. (j83p may be rewritten as 



5Dii^{z) ^ 



D. 



off 



(84) 



where SDib{z) stands for the quasi-one-dimensional 
weak localization correction, and ^id = nvaDcft is the 
exact localization lengt h^^i^i . Eq. agrees with the 
leading z/^iD-expansion of the local diffusion coefficient 
given in Ref. It thereby justifies that for z > a the 
medium bar displays the quasi-one-dimensional (inter- 
face) weak localization. Indeed, at early times t < Dq ja} 
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incident photons explore a region of size c? neighboring 
to the interface, approaching a uniform distribution in 
the transverse direction. At later times they diffuse as 
in a quasi-one-dimensional medium. Technically, start- 
ing from quasi-one-dimensional a model by performing 
the one- loop calculation one finds Eq. (I84|) . Importantly, 
from Eqs. ([5^ and (IM)) we find that within the boundary 
layer 2: < ^ weak rather than strong localization occurs 
even in the static limit. 



C. CBS line shape: crossover from 2D weak to 
quasi- ID strong localization 

In this part we turn to investigate effects of local diffu- 
sion on the CBS line shape. We will consider a medium 
illuminated by light of frequency fl parallel to the bar, 
and calculate the angular resolution of the backscattered 
light intensity a{6) near the inverse incident direction. 
Since the bar is wide enough so that I ^ a I ) 
a large parametric region: X/a < 9 < X/l is opened. 
Below we pay particular attention to the line shape at 
< ^ A//. (Notice that the line shape is symmetric 
with respect to 6 — .) 

It is well known that the backscattered light intensity 
may be decomposed into the background ao and the co- 
herent part ac{0) according to 



a{e) = aQ + aa{9) . 



(85) 



Here 



drdr' 1 



-3^°(r,r') 



(86) 



cos[qi-(r-r')]3^^(r,r') 



where q± = 217 sin(6'/2) w 2tt6/\ because of ^ 1 , and 
the overall normalization factor is omitted. First of all, 
it is easy to show that ao = ac(0) and therefore only 
the coherent part ac{9) which determines the line shape 
will be studied below. Inserting the Fourier transform, 
namely Eq. (|76p into ac{d) we arrive at 



dzdz' e-^y'? {z,z') 



/o Jo 



(87) 



where the propagator y^^{z, z') solves Eq. (|77)) . 



1. Signatures of 2D weak localization 

The interfering optical paths penetrate into the 
medium of a depth ~ 1^ ■ Let us first study the 
CBS line shape in the region: n/a ^ q±^ < l^^ . Be- 
cause of the condition: / < gj^ <C a the CBS line 
shape is mainly responsible for by photons which dif- 
fuse around the interface, i.e., Z z <C a and thereby 



undergo two-dimensional weak localization. Indeed, the 
first term of Eq. ([5T|) is much smaller due to the condition 
\q^a\/'K ^ 1 . Setting the ultraviolet cutoff N ^ a/{nl) 
(as Eq. (pn)) ) we may approximate Eq. ([5T|) by 



SD{z) 

Dn 



a/{TTl) 

a X ^ 

aqj_/7T 



With the substitution of such weak localization correc- 
tion into Eq. ((77)) we find 



^0 



1- — Inlg^^l 



n/a<$:qj_<l ^ 



(89) 



According to Eq. (|89p the conventional triangular peak 
described by the factor l — 2lq± is enhanced by a logarith- 
mic factor (dashed line in Fig.[T]). Eq. ([88]) indicates that 
in the region: ir/a q± ^ the local diffusion is of 
minor importance. It is the two-dimensional bulk weak 
localization that is responsible for such a logarithmic en- 
hancement. Thus in such region the scaling theory, still, 
is applicable. 



2. Signatures of quasi-lD strong localization 

At the exact backscattering direction, i.e., = in- 
terfering optical paths penetrate into the medium bulk: 
z ^ £, where quasi-one-dimensional (bulk) strong local- 
ization states are formed. In contrast to ir/a <^ q± ^ 
for q± ^ the local diffusion plays crucial roles and 
thus strongly affects the CBS line shape. Indeed, de- 
spite of the nonperturbative nature of strong localization 
the backscattering light intensity at = may be easily 
found provided that in the region z ^ ^ the local diffu- 
sion equation is still valid^^. For q± = from Eq. ([77)1 
one may find 



dz — -— 

D{z) 



I 



D. 



off 



(90) 



where in the second equality the substitution of Eq. (183 
is made. This immediately gives 



ac(0) 



/3 



(91) 



Eq. (PT|) shows that although due to two-dimensional 
bulk weak localization the CBS line shape develops a log- 
arithmic singularity at n/a < q± ^ , the singularity 
is cut off at gj^ = (Fig. [1]) where quasi-one-dimensional 
strong localization occurs in the bulk. 

It is in order to remark that formally there is a region 
q±^ < 1 where the rounding due to the local diffusion oc- 
curs (Yet, the detailed rounding form depends on D(z) in 
the region z ^ ^ , to find which is beyond the present per- 
turbative treatise.), however, it is unobservable because 
the finite bar width renders <IC 7r/a . 
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FIG. 1: (Color online) Coherent backscattering intensity 
(solid circle) in unit of l^/Do versus q± = niv/a for the 
medium bar with a/l = 100 . Upper panel: The logarith- 
mic enhancement-dashed line-is cut off at = . From top 
to bottom the parameter Z/A is 2 (red), 10 (green), 20 (purple) 
and (x> (blue). Lower panel: The magnification of the upper 
panel for lq± closed to (The blue curve is not plotted.). 



Finally we anticipate that the predicted CBS line shape 
is qualitatively correct for I > X , though the analytical 
result here is obtained for Z ^ A . 



VI. CONCLUSIONS 

For light propagation in fully disordered media a super- 
symmetric field theory is presented. The supersymmetric 
a model described by Eq. (fT5)) may be applied to bulk 
(infinite) media for studies of optical localization transi- 
tion. In this direction it may serve as an alternative to the 
replica field theoryi. However, the supersymmetric field- 
theoretic formalism turns out to be far more powerful as 
propagation of incident light in semi-infinite media— con- 
cerned, which is the subject of this paper and is closely 
related to the coherent backscattering phenomenon. 

Differing from infinite medium in the presence of the 
vacuum-medium interface C (The interface may bear ar- 
bitrary geometry but must be smooth over the scale of 
the mean free path.) the supermatrix field (locally) sat- 
isfies the radiative boundary condition Eq. (PD|) . Ac- 



cordingly, the bare diffusion constant acquires a position- 
dependent wave interference correction namely 



D{r;Lu) = DQ + 6D{r;u;). 



(92) 



which roots in the incomplete constructive interference 
(weak localization) near the interface. Thus, we jus- 
tify the (static) local diffusion equation, originally pro- 
posed in Ref. |1^ as well as its dynamic (i.e., uj ^ 0) 
generalizationii. Most importantly, for (almost) trans- 
parent interface, i.e., To{r) « l,r G C, the weak lo- 
calization correction SD(r]uj) vanishes at the interface. 
This immediately shows that the radiative boundary con- 
dition is protected against wave interference effects, and 
constitutes an explicit proof that no scaling hypothesis 
might exist in the extrapolation layer of thickness ~ / . 
Therefore, the present work supports the criticism of 
Refs. [HB 

on the earlier theoretical proposal^. 
In the present work the static limit of the wave interfer- 
ence (weak localization) correction namely 5D{r ; ijj ^ 0) 
is explicitly calculated for the two-dimensional semi- 
infinite medium with a finite width a (the bar geom- 
etry), where 6D{r;uj 0) solely depends on the dis- 
tance from the interface z . For Z <C a ^ le"^ " a 
dimensional crossover of the wave interference correction 
5D{z) = SD{r;uj 0) is found. Indeed, SD(z) displays 
two-dimensional weak localization at z <C a with a loga- 
rithmic dependence on z , while displays one-dimensional 
weak localization at a <C z (^ ^) . Furthermore, for the 
latter region it is not difficult to generalize Eq. ([M]) to 
higher order loop corrections which reads out as 



SDiuiz) 



D. 



off 



Z 

1^ 



n=2 



Z 



(93) 



This-at the perturbative level-formally confirms the re- 
sult of Ref. llj for one-dimensional geometry. Notice that 
the unimportant numerical expansion coefficients c„ may 
vary depending on the strict/quasi- one-dimensional ge- 
ometry. 

For wider medium bar such that a ^ le'^ (e-g-j 
infinite medium plane) Eq. ([82]), in fact, is still applicable 
except that the one-dimensional contribution vanishes. 
That is, 



5D{z) 
Do 



1 , £ 
n^i^Do I 



(94) 



This suggests that in the medium there exists a bound- 
ary layer of thickness le^ outside which two- 
dimensional strong localization occurs. Surprisingly, in- 
side the layer the diffusion coefficient logarithmically de- 
pends on the distance from the interface. How to repro- 
duce this logarithmic dependence by the self-consistent 
diagrammatical niethod^'^ is unclear. 

Finally, it should be stressed that the present field- 
theoretic justification of local diffusion is perturbative. 
Therefore, the validity of such a concept in the nonper- 
turbative strong localization region remains an impor- 
tant question. This problem is far beyond the scope of 
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this paper and will be addressed in a forthcoming pa- 
per, especially the issue how Eq. ([M)) is extended to the 
nonperturbative region: z ^ ^ . It is well known that in 
Faraday-active medium the one-loop weak localization 
may be strongly suppressed'^'^ . Therefore, to take into 
account such medium within the present field-theoretic 
formalism remains another important problem. It is also 
interesting to generalize the present field-theoretic for- 
malism to include medium gain'^'*. These issues are left 
for future work. 



Acknowledgments 

I am deeply grateful to A. Altland, T. Micklitz for use- 
ful conversations, especially M. R. Zirnbauer for several 
fruitful discussions. Work supported by Transregio SFB 
12 of the Deutsche Forschungsgemeinschaft. 



where (7+ (C~ ) stands for the curve infinitesimally closed 
to C from the V+ (V-) side. Taking Eq. (|A5p into account 
we may rewrite Eq. ([3]) as 

[n\-H] GPA^y)- [ dv"B{vy)GfA^",v') 

J C 

(A7) 



= 0, reC, 



which is supplemented by the boundary condition 
9„(r)G« (r,r')|rec = 0. 

i?(r,r') (r,r' G C) consists of the real and imaginary 
part. The former is small and may be absorbed into fl^ 
renormalizing e(r) . It is thus ignored. In contrast, the 
latter is important determining the analytical structure. 
Taking this into account we prove Eq. ([^0]) for the re- 
tarded (and similarly for the advanced) Green function. 



APPENDIX A: PROOF OF EQ. ^ 



APPENDIX B: SIMPLIFICATION OF THE 
COUPLING ACTION Filter [Q] 



Assuming that r , r" G V_ and r' e V+ , from the 
Helmholtz equation ([T]) we obtain: 



G^. (r, r'){V^ + n^il + 6(r)]}[g^. (r, r")]' 



(5(r-r")G« (r,r') (Al) 



and 



(r, r")]* {V^ + nl[l + 6(r)]}G«. (r, r') = . (A2) 

Subtracting Eq. (jA2p from Eq. (IA1[) gives 

V . {[Vg,t(r,r")]G«.(r,r') - [5jt.*(r, r")]VG«. (r, r')} 
= 5(r-r")G^.(r,r'), (A3) 

where V acts only on r . Noticing that g^2 (r, i"")|rec = 
and [g^2(r, r")]* = g^2(r",r), with r e V_ integrated 
out we find 

G«.(r",r')= / dr9„(r){5^.(r",r)}G^2(r,r'). (A4) 
Jc 

Here (?„(r) stands for the normal derivative at r with n(r) 
pointing to V-j- . Taking the derivative we obtain: 

an(r")G^^(r",r')= / dri?(r",r)G«.(r,r') (A5) 
Jc 



with B(r" , r) following the definition of Eq. ((2T|) . 

Now suppose that r is shuffled to C from the V_ side. 
Let us integrate out Eq. ([3|) over the line element along 
an infinitesimal piece of a curve passing from V_ to V+ . 
In doing so we obtain: 

9„(r)Gg.(r,r')|rec+ -a„(r)Gg.(r,r')|rec- =0, (A6) 



For the moment let us suppress the indices i and k± , 
and decompose Q according to 



[Q±,A]^0, {Q||,A} = 0. 



(Bl) 



Taking such decomposition into account we obtain: 

strln(l + aAQ) (B2) 



strln (l + aAQ\\) + strln ( 1 + j 



1 



aAQi 



■ aAQj 



Upon Taylor expanding the second logarithm only the 
even order terms contribute. Thus, Eq. (jB2p may be 
rewritten as 

strln(l + aAQ) 

= strln (1 + aAQii) + i strln (1 + aAQj_ 

2 \ l + aAQji 

i strln 1 1 - ^ . ^ aAQj 

= i strln {(l + aAO[|)' + (aQi)'} 
= i strln (2 + 2q;^ + 4q!A(3||) 

= istrln(2-ro + ToAQ||) , (B3) 

where in deriving the third equality we use the identity 
Q| + Ql = 1 J E^nd in deriving the last two equalities we 
use the identity str In 1 = . 

Restoring the index i and substituting Eq. (|B3p into 
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Eq. ([291) we obtain 

FinteAQ] = ^strln[2-roW+ToWAQ,||] 

i 
i 

-J^^ToWstr (Ag,) , 



(B4) 



where in the second hne we take advantage of strong 
coupHng, i.e., Ql 1 , and in the third cquahty we use 
the identity str {AQj_) = . 



APPENDIX C: THE CHARGE-CONJUGATION 
SYMMETRY OF W 



giving kB'f = B and B'^k = B . The first relation may 
be rewritten as 

B*k = CoBCj. (C2) 

Substituting the second relation into Eq. (jC2p we obtain: 

B*k = CoS*^fcCj (C3) 

giving Bk = CoB^kC^ . Noticing the relation: 
Cq/cCJ = k one finds 



B — Cq_B^C(F = B 



(C4) 



and thus justifies the charge-conjugation symmetry of 
W . 



The charge-conjugation symmetry is irrespective of 
fast-slow mode separation and therefore we ignore the 
subscript > (<) . Substituting the parametrization of W 
into Eq. (j48p we obtain: 



B 
B 



B^kV 
kB^ ) 



(CI) 



APPENDIX D: PRESERVATION OF WARD 
IDENTITY 

Using integral by parts we transform Eq. (j65p into 
(noticing that dpW{r)\p^oora = 0) 



Fi[W] 



oo pa 



'0 JO 

Using the contraction rule we obtain: 



dz / dp ^^2Do str d^WW^ + {dW wf + [dwf W'^ 



iLO str W 



oc ca 



TTiyDo / dzi dpi {stT[k{l + A){l-T3)W{r)k{l^A){l + T3)W{r')]sti[d^W{ri)W^{ri)]) 
Jo Jo 



F2 



(Dl) 



(D2) 



OC pa 



la + TTiyDo / dzi / dpi (str [kil + A)(l - T3)W{r)k{l - A)(l + T3)W{r')] str [d^ W {ri)W{ri)W (ri)VK(ri)] ) , 



/O ^0 

where the overbrace fixes the contraction and 

^•00 pa 



la = TTuDa j dzijdpi\{sirk{l + K){l-T:i)W{T)k{l-A){l + T3)W{T')sivd''W[ri)W\ri))p^ + 



(str/c(l + A)(l - T3) W{Y)k{l ~ A)(l -I- T3)W{r') aird^W {y{)W^{yi)) 



F2 



= Do J dzij dpi}^dlV{r',ri;Lj){stT[k{l + A){l-T3)W{r)k{l-A){l+T3)W\ri)])j,^ + 
d^Vir, n; w)(str [k^l A)(l - T3)W^iri)kil - A)(l + T3)M^(r')])^^ 

Likewise, we also obtain: 

p 00 pa 

h = ^ y^dziy^dpi (str[fc(l + A)(l-T3)M^(r)fc(l- A)(l + T3)W^(r')]str[T4^4(ri)])^^ 

poo pa 1' 

= j dziydpi|l?(r',ri;cc;)(str[fc(l + A)(l-r3)M^(r)fc(l- A)(l + T3)W^3(ri)])^^ + 
2?(r, Vi-Lo) (str [fc(l + A)(l - T3)W\vi)k{l - A)(l + T3)W{v')])^ 



(D3) 



(D4) 
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Noticing Eq. ((67|) we find that la + h exactly cancels the 
first two terms of Eq. . 

APPENDIX E: DERIVATION OF EQ. (f82|) 

Let us introduce the function: f{z) = X]n>i ^ 
g-2n7rz/aj^^j^^-j . Taking its derivative we obtain: 

/'(z) = -Y: e-'-'^^ = - / , . (El) 

a ^ — ' a 1 — e ^^^i'^ 



On the other hand, the low-energy diffusion occurs on the 
scale ^ I . Over this scale the interface where /(z) van- 
ishes is smeared. Therefore, without loss of any physics 
we may reformulate the boundary condition as /(/) = . 
Taking it into account and integrating out Eq. (|E1[) we 
obtain 7r/(z) = ln{(l-e-2^^/'')/(l-e-2'^'/'')} for z > / , 
which gives J(z) k, tt^^ ln(z/Z) for I ^ z a justifying 
Eq. dMl). 
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